Domain wall in a chiral p-wave superconductor: a pathway for electrical current 
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Superconductors with p x ± ip y pairing symmetry are characterized by chiral edge states, but these 
are difficult to detect in equilibrium since the resulting magnetic field is screened by the Meissner 
effect. Nonequilibrium detection is hindered by the fact that the edge excitations are unpaired 
Majorana fermions, which cannot transport charge near the Fermi level. Here we show that the 
boundary between p x + ip y and p x — ip y domains forms a one-way channel for electrical charge. We 
derive a product rule for the domain wall conductance, which allows to cancel the effect of a tunnel 
barrier between metal electrodes and superconductor and provides a unique signature of topological 
superconductors in the chiral p-wave symmetry class. 

PACS numbers: 74.20.Rp, 74.25.fc, 74.45. +c, 74.78.Na 



Chiral edge states are gapless excitations at the bound- 
ary of a two-dimensional system that can propagate in 
only a single direction. They appear prominently in the 
quantum Hall effect [TJ |2] : The absence of backscatter- 
ing in a chiral edge state explains the robustness of the 
quantization of the Hall conductance against disorder. 
Analogous phenomena in a superconductor with broken 
time reversal symmetry are known as the spin quantum 
Hall effect [31 IH |S] and the thermal quantum Hall effect 
[5J [7] , in reference to the transport of spin and heat along 
chiral edge states. 

Unlike the original (electrical) quantum Hall effect, 
both these superconducting analogues have eluded obser- 
vation, which is understandable since it is so much more 
difficult to measure spin and heat transport than electri- 
cal transport. Proposals to detect chiral edge states in a 
superconductor through their equilibrium magnetization 
are hindered by screening currents in the bulk, which 
cancel the magnetic field (Meissner effect) [5J [HI EH EQ • 

Here we show that the boundary between domains of 
opposite chirality (p x ± ip y ) in a chiral p-wave supercon- 
ductor forms a one-way channel for electrical charge, in 
much the same way as edge states in the quantum Hall 
effect. This is not an immediate consequence of chiral- 
ity: Since the charge of excitations in a superconductor 
is only conserved modulo the Cooper pair charge of 2e, 
the absence of backscattering in a superconducting chiral 
edge state does not imply conservation of the electrical 
current. Indeed, one chiral edge state within a single 
domain has zero conductance due to electron-hole sym- 
metry. We calculate the conductance of the domain wall, 
measured between a pair of metal contacts at the two 
ends (see Fig. [IJ, and find that it is nonzero, regardless 
of the separation of the contacts. 

Our analysis is generally applicable to so-called class-D 
topological superconductors [HUH], characterized by the 
presence of electron-hole symmetry and the absence of 
both time-reversal and spin-rotation symmetry. It can be 
applied to the various realizations of chiral p-wave super- 
conductors proposed in the literature (strontium ruthen- 
ate [11], superfluids of fermionic cold atoms [HI [15], and 
ferromagnet-superconductor heterostructures |16l 117] ) . 




FIG. 1: Superconducting strip divided by a domain wall 
(dashed line, length W) into domains with p x ± ip y symme- 
try. The edge states ^l, of opposite chirality in the two 
domains are indicated by red arrows. These unpaired Majo- 
rana modes can carry heat current between contacts Nl and 
Nr, but no electrical current. A normal metal electrode Ni 
at voltage Vi injects charge into the domain wall, which is 
detected as an electrical current I2 at the other end N2. In 
an alternative measurement configuration (indicated in blue) , 
contact N2 measures a voltage V2 without drawing a current. 



We start from the Bogoliubov-De Gennes equation, 



H — Ep 
At 



-m + E f 



E 



(1) 



for coupled electron and hole excitations u(r),v(r) at en- 
ergy E above the Fermi level Ep. The single-particle 
Hamiltonian is Hq = (p + eA) 2 /2m + U, with p = 
—ihd/dr the momentum, A(r) the vector potential, and 
U(r) the electrostatic potential. The dynamics is two- 
dimensional, so r = (x, y), p = (PxiPy)- The pair poten- 
tial A has the spin-polarized-triplet p-wave form [18 



A = (2p F ) I (v-P + P-V), 



(2) 
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in terms of a two-component order parameter r) — 
iVxtVy)- The two chiralities p x ± ip y correspond to 
T)± = A e J *(l,±i), with A the excitation gap and 
the superconducting phase. Since = —A*, a solution 
(u,v) ofEq. Q at energy 75 is related to another solution 
(v*,u*) at energy — E (electron-hole symmetry). A do- 
main wall along x — 0, with a phase difference <f> between 
the domains, has order parameter |19l I20j 

Vx (x) = A Q [e^ 2 cosx(x) + e l ^ 2 sm X (x)}, (3a) 
rj y {x) = iA [e- l4,/2 cos X (x) - e 1 ^ 2 sm X (x)}, (3b) 

The function x( x ) increases from to ir/2 over a coher- 
ence length £o — ^vf/Aq around x = 0. 

At energies E below Ao the excitations are nondcgcn- 
erate chiral edge states and ^r circulating in opposite 
directions in the two domains [5J [5TJ [551 HHJ- (See Fig. 
[I]) At the domain wall the two states mix, so that an ex- 
citation entering the domain wall in the state *i>™ or 
can exit in either of the two states \E , ™ t and ^ R ut . We 
first analyze this edge state scattering problem between 
contacts Nl and TVr, and then introduce the contacts TVi 
and 7V 2 to the domain wall. 

The edge state excitations have fermionic annihila- 
tion operators j(E) = l{E) ^ r{E)~) , which satisfy the 
electron-hole symmetry relation 

1 (E) = 1 \~E). (4) 

At zero energy one has 7 = ~f , so these are Majo- 
rana fermions [18j . The unitary scattering matrix S(E) 
relates incoming and outgoing operators, 7° ut (7?) = 
S(E)~f ln (E). Electron-hole symmetry for both 7 111 
and 7 out requires S(E)y n (E) = -y in (E) S^ (- E) , hence 
S(E) = S*(—E). The zero-energy scattering matrix 
S(0) = Sdw of the domain wall is therefore a real uni- 
tary, or orthogonal, matrix. We may parametrize it by 

adw ~ \(-l)P +1 smiP (-l)PcosVy z ' [ > 

in terms of a mixing angle tp and a parity index p£ {0 , 1 } . 

The mixing angle ip — k y W is determined by the phase 
accumulated by the pair of chiral Majorana modes, as 
they propagate with wave number ±k y along the domain 
wall of length W. The dispersion relation E(k y ) of the 
Majorana modes was calculated in Ref. [22 , for a step 
function order parameter at x — 0, including also the 
effect of a tunnel barrier U = UqS(x) (tunnel probability 
D, zero magnetic field). By equating E(k y ) — and 
solving for k y we obtain the mixing angle 

ip = k F wVDcos((j)/2). (6) 

The mixing angle can in principle be measured through 
thermal transport between contacts Nl and Nr, since 
the heat current through the domain wall is cx sin 2 ip. 
In what follows we consider instead a purely electrical 
measurement of transport along the domain wall, that 



(as we shall see) is independent of the degree of mixing 
of the Majorana modes. 

The measurement that we propose consists of the injec- 
tion of electrons from contact N\ at voltage V\ (relative 
to the superconductor) and the detection at contact -/V 2 . 
We consider two detection schemes: In the first scheme 
contact N2 is kept at the same potential as the supercon- 
ductor and measures a current I2, leading to the nonlocal 
conductance Gi 2 = h/Vi- I n the second scheme contact 
N2 is a voltage probe drawing no net current and mea- 
suring a voltage V%. The ratio R12 — V2/I1, with I\ the 
current entering the superconductor through contact N\ , 
is the nonlocal resistance. The two nonlocal quantities 
are related by R12 — G12/G1G2, with d = \h/Vi\ the 
contact conductance of electrode Ni (measured with the 
other contact grounded). 

We take the zero-temperature and zero-voltage limit, 
so that we can use the zero-energy scattering matrix to 
calculate the various conductances. The scattering prob- 
lem at contact Ni involves, in addition to the Majorana 
operators 7 = (7l,7h), the electron and hole annihila- 
tion operators a n and b n in mode n = 1, 2, ... TV. These 
are related by b n (E) — a^—E). The even and odd com- 
binations 7„ , defined by 




satisfy the same electron-hole symmetry relation Q as 
1l,1r, and therefore represent Majorana fermions at 
E = 0. We denote j n = (7^,7^) and collect these op- 
erators in the vector T = (71, 72, • ■ • In)- The scattering 
matrix S± of contact N\ relates incoming and outgoing 
operators, 

(?)_"*(?), <»> 

Electron-hole symmetry implies that Si is (27V + 2) x 
(2N + 2) orthogonal matrix at zero energy. Similarly, 
the zero-energy scattering matrix S2 of contact iV 2 is a 
(2N' + 2) x (2N' + 2) orthogonal matrix. (The number 
of modes is N, N' in contacts N\, N2 respectively.) 
The 2N' x 27V transmission matrix 

tai = t' 2 S dw ti = t^e^h (9) 

from contact Ni to N 2 is the product of the 2 x 27V 
submatrix t\ of S\ (transmission from N\ to the domain 
wall) , the 2x2 scattering matrix Sdw (transmission along 
the domain wall), and the 27V' x 2 submatrix t' 2 of 5 2 
(transmission from the domain wall to 7V 2 ). 

The total transmission probability T eei summed over 
all modes, of an electron at contact TVi to an electron at 
contact 7V 2 is given by 

T ee = \YxUH\ Y U{\ + ^ z )Uh 21 U{\ + S,) (10) 
= iTr4(l-S y )t 21 (l-S I/ ), (11) 
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where we have defined the direct sums U = u © u ■ ■ ■ © u, 
Si = <7j © <Ji ■ ■ ■ © <ii and we have used that ua z u^ = 
—<T y . Similarly, the total electron-to-hole transmission 
probability The reads 

T fee = iTr4(l + S s/ )t 21 (l-Sj / ). (12) 

Since 1% = (e 2 / h)Vi(T ee — The), the nonlocal conduc- 
tance takes the form 



G 12 = {e 2 /h)^Trt^E y t 21 E y 



(13) 



We have used that t 21 = iji an d Tri^Xj^i = (being 
the trace of an antisymmetric matrix) . The nonlocal re- 
sistance can be written in a similar form upon division 
by the contact conductances, 

Ri2 = -£yr, G i = ( e 2 /M^Tr(l-E„rfs,r.:). (14) 

We will henceforth set e 2 /h to unity in most equations. 

Substitution of Eq. ^ into Eq. (13 1 gives the conduc- 
tance 



G12 — \^T\S^ w T 2 Sdw, 



(15) 



in terms of the 2x2 matrices 7{ — t{S y ti, T 2 = t' 2 T, y t' 2 . 
We now use the identity 

TtAxA 2 = \ (Tr Aya y ) {Tr A 2 a y ) , (16) 

valid for any pair of 2 x 2 antisymmetric matrices A\, A 2 . 
Taking A x = Ti, A 2 = Sj w T 2 S dw we arrive at 

G12 = (-l) p aiQ!2, a, = yYx%o y , (17a) 
R12 = (-1)W 2 , p i = on/G i , (17b) 

since Tr Sj , n T 2 Sd w (J v = {-l) pr TT:T 2 u y in view of Eq. 

Eq. ( |l~7| ) expresses the nonlocal conductance and resis- 
tance in terms of the scattering matrices S±, S 2 of the two 
contacts Ni,N 2 . The scattering matrix Sdw of the do- 
main wall enters only through the parity index p, and not 
through the mixing angle if). That the transferred charge 
depends only on a parity index is a generic feature of a 
single-mode scattering problem with class D symmetry 
El E3 H3 EH |28]. Quite generally, p counts the num- 
ber (modulo 2) of zero-energy bound states, which in our 
case would be trapped in vortices in the domain wall. 

A measurement of the domain wall conductance would 
have several characteristic features: Most prominently, 
the conductance is zero unless both contacts N\ and N 2 
are at the domain wall; if at least one contact is moved 
away from the domain wall, the conductance vanishes 
because a single Majorana edge mode cannot carry an 
electrical current at the Fermi level [29 . This feature 
would distinguish chiral p-wave superconductors (sym- 
metry class D) from chiral <i-wave superconductors (sym- 
metry class C), where the Majorana edge modes come in 
pairs and can carry a current. The chirality itself can 
be detected by interchanging the injecting and detecting 
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FIG. 2: Solid curves: probability distributions of the nonlo- 
cal conductance G12 (in units of e 2 /h) and nonlocal resistance 
R12 (in units of h/e 2 ). These are results for a random distri- 
bution of the 4x4 orthogonal scattering matrices 5i and S2 ■ 
The dashed curve shows the narrowing effect on P(Gi2) of a 
tunnel barrier in both contacts (tunnel probability r = 0.1). 
In contrast, P{R\2) is not affected by a tunnel barrier. 



contacts: only one choice can give a nonzero conduc- 
tance. While vortices trapped in the domain wall can 
change the sign of the conductance (through the parity 
index p) , other properties of the domain wall have no ef- 
fect on G\ 2 . In particular, there is no dependence on the 
length W. 

To illustrate these features in a model calculation, we 
consider the case of two single-mode contacts (N = N' = 
1) coupled to the domain wall through a disordered in- 
terface. We model the effect of disorder using random 
contact scattering matrices S\ and S 2l drawn indepen- 
dently with a uniform distribution from the ensemble of 
4x4 orthogonal matrices. In the context of random- 
matrix theory [30] . uniformly distributed ensembles of 
unitary matrices are called "circular" , so our ensemble 
could be called the "circular real ensemble" (CRE) — to 
distinguish it from the usual circular unitary ensemble 
(CUE) of complex unitary matrices [31] . 

Using the expression for the uniform measure on the 
orthogonal group E3 [33] , we obtain the distributions of 
the parameters ai and (3i characterizing contact Nf 



P(a) = l-\a\, P(/3) = (l + 



\a\,\P\<l. (18) 



The distribution of the nonlocal conductance G\ 2 — 
( — l) p aiO!2, plotted in Fig. [2j then follows from 

P(G 12 ) = J d ai J da2<5(Gi2-aia2)P(ai)P(a 2 ) 

= 4|G 12 |-4-2(l + |G 12 |)ln|G 12 |, |G 12 | < 1. 

(19) 

(There is no dependence on the parity index p because 
P is symmetric around zero.) The distribution of the 
nonlocal resistance R\ 2 — (~l) p (3iP 2 follows similarly 
and as we can see in Fig. [2] it lies close to P{G\ 2 ). 

The difference between the two quantities Gi 2 and R\ 2 
becomes important if the contacts between the metal and 
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the superconductor contain a tunnel barrier. A tunnel 
barrier suppresses G12 but has no effect on R\2- More 
precisely [33], any series resistance in the single-mode 
contacts N\ and N2 which does not couple electrons and 
holes drops out of the nonlocal resistance R\2- This re- 
markable fact is again a consequence of the product rule 
(16 1, which allows to factor a series conductance into a 
product of conductances. A tunnel barrier in contact i 
then appears as a multiplicative factor in a* and Gi , and 
thus drops out of the ratio ft = a*/ Gi determining R 12 ■ 
To demonstrate the effect of a tunnel barrier (tunnel 
probability t), we have calculated the distribution of a 
using the Poisson kernel of the CRE [52] , with the result 



P(a,r) = 



[r+ (1 -T)\a\} 3 [t + (1-t> 



212 ' 



(20) 



The distribution of f3 remains given by Eq. ( 18 1, indepen- 
dent of r. The dashed curves in Fig. [2] show how the re- 
sulting distribution of the nonlocal conductance becomes 
narrowly peaked around zero for small r, in contrast to 
the distribution of the nonlocal resistance. 

Among the various candidate systems for chiral p- 
wave superconductivity, the recent proposal [T^] based on 



the proximity effect in a semiconducting two-dimensional 
electron gas seems particularly promising for our pur- 
pose. Split-gate quantum point contacts (fabricated with 
well-established technology) could serve as single-mode 
injector and detector of electrical current. The chiral- 
ity of the superconducting domains is determined by the 
polarity of an insulating magnetic substrate, so the loca- 
tion of the domain wall could be manipulated magneti- 
cally. The appearance of a nonlocal signal between the 
two point contacts would detect the domain wall and the 
disappearance upon interchange of injector and detector 
would demonstrate the chirality. 

As a direction for further research, we note that do- 
mains of opposite chirality are formed spontaneously in 
disordered samples. Since, as we have shown here, do- 
main walls may carry electric current, a network of do- 
main walls contributes to the conductivity and may well 
play a role in the anomalous (parity violating) current- 
voltage characteristic reported recently [31] . 
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and by an ERC Advanced Investigator Grant. 



1] B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 
2] M. Buttiker, Phys. Rev. B 38, 9375 (1988). 
3] G. E. Volovik and V. M. Yakovenko, J. Phys. Cond. Mat- 
ter 1, 5263 (1989). 
1] T. Senthil, J. B. Marston, and M. P. A. Fisher, Phys. 

Rev. B 60, 4245 (1999). 
5] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
3] T. Senthil and M. P. A. Fisher, Phys. Rev. B 61, 9690 
(2000). 

7] A. Vishwanath, Phys. Rev. Lett. 87, 217004 (2001). 
3] M. Matsumoto and M. Sigrist, J. Phys. Soc. Japan 68, 
994 (1999). 

9] H. J. Kwon, V. M. Yakovenko, and K. Sengupta, Synth. 

Met. 133, 27 (2003). 
3] J. R. Kirtley, C. Kallin, C. W. Hicks, E.-A. Kim, Y. Liu, 
K. A. Moler, Y. Maeno, and K. D. Nelson, Phys. Rev. B 
76, 014526 (2007). 
1] C. Kallin and A. J. Berlinsky, J. Phys. Cond. Matter 21, 
164210 (2009). 

2] A. Altland and M. R. Zirnbauer, Phys. Rev. B 55, 1142 
(1997). 

3] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- 

wig, Phys. Rev. B 78, 195125 (2008). 
1] S. Tewari, S. Das Sarma, C. Nayak, C. Zhang, and P. 

Zoller, Phys. Rev. Lett. 98, 010506 (2007). 
5] M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 

103, 020401 (2009). 
3] J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 

larXiv:09072239l 

[17] P. A. Lee, |arXiv:0907.2681| 

[18] A. Stern, F. von Oppen, and E. Mariani, Phys. Rev. B 

70, 205338 (2004). 
[19] M. Sigrist and D. F. Agterberg, Prog. Theor. Phys. 102, 



965 (1999). 

[20] A. Bouhon and M. Sigrist, [arXiv:0909.3535| 

[21] Yu. S. Barash, A. M. Bobkov, and M. Fogelstrom, Phys. 

Rev. B 64, 214503 (2001). 
[22] H. J. Kwon, K. Sengupta, and V. M. Yakovenko, Eur. 

Phys. J. B 37, 349 (2004). 
[23] M. Stone and R. Roy, Phys. Rev. B 69, 184511 (2004). 
[24] B. Beri, J. N. Kupferschmidt, C. W. J. Beenakker, and 
P. W. Brouwer, Phys. Rev. B 79, 024517 (2009). 
L. Fu and C. L. Kane, Phys. Rev. Lett. 102, 216403 
(2009). 

[26] A. R. Akhmerov, J. Nilsson, and C. W. J. Beenakker, 

Phys. Rev. Lett. 102, 216404 (2009). 
[27] B. Beri, Phys. Rev. B 79, 245315 (2009). 



[25] 



[28] K. T. Law, P. A. Lee, and T. K. Ng, arXiv:0907.1909 



[29] That the nonlocal conductance vanishes if one of the two 
contacts couples only to a single domain, can be seen 
directly from Eq. (17 1: If, say, contact 1 couples only to 



the right domain, then only the 2, 2 element of T\ can be 
nonzero, but since this matrix is antisymmetric the 2, 2 
element must also vanish and 71 must be zero identically. 
This implies ot\ = 0, hence G12 = 0. 
[30] M. L. Mehta, Random Matrices (Elsevier, Amsterdam, 
2004). 

[31] The name "circular orthogonal ensemble" (COE) might 
be more appropriate for the ensemble of uniformly dis- 
tributed orthogonal matrices, but this name is already in 
use for the ensemble of unitary symmetric matrices. 

[32] B. Beri, Phys. Rev. B 79, 214506 (2009). 

[33] For a detailed calculation we refer to the appendices. 

[34] H. Nobukane, A. Tokuno, T. Matsuyama, and S. Tanda, 
larXiv:0906. 36441 



5 



APPENDIX A: AVERAGES OVER THE 
CIRCULAR REAL ENSEMBLE 



To calculate the distributions (18) of the parameters 



cti and Pi we need the probability distribution of the 4x4 
scattering matrix Si of contact i = 1,2 in the CRE. We 
may either work in the basis of electron and hole states, 
as in Rcf. [27 , or in the basis of Majorana states. Here 
we give a derivation of Eq. (18 1 using the latter basis 



(which is the basis we used in the main text). 

A 4 x 4 orthogonal scattering matrix has the polar 
decomposition 



S = 



S 



c 



) \(-l) p+1 C (~l) p S 



c 



( cos ipi 

I COS tp2 



s 



( sin tpi 
I sin ip2 



(A2) 



in terms of six real angles. We need the uniform mea- 
sure on the orthogonal group, which defines the proba- 
bility distribution in the circular real ensemble (CRE). 
This calculation proceeds along the same lines as in Ref. 
|27j (where a different parametrization, in the electron- 
hole basis, was used). The result is that the angles 
</>i, </>2, 03, 04 are uniformly distributed in (0,27r), while 
the angles ipi , ip2 have the distribution 

P(ipi,ip 2 ) = jlcos^i-cos 2 ^!, 0< Vi>2 < 7T. (A3) 

We can now obtain the joint distribution P(a.i,Gi) of 
the injection (or detection) efficiency and the (dimen- 
sionless) contact conductance Gi of contact i. (We drop 
the label i for ease of notation.) By definition, 



^Tr t(Jyt T Oy 



a = % Ltt<T v t' o y = cos ipi cos 1p2, 
G = 1 — ^Tr rUyT 1 a y = 1 — sin ip\ sin ip2- 
Notice the trigonometric inequality 

< \a\ < G < 2 - \a\. 



(A4) 
(A5) 

(A6) 

By averaging over the CRE we find, remarkably enough, 
that the joint distribution of a and G is uniform when 
constrained by this inequality, 

fir />7r 

P(a,G)= / #! / dfaPtyuih) 
Jo Jo 

x 6(a — cos ■01 cos'02)^(G — 1 + sin-0i sin 1P2) 

' 1/2 if < |a| < G <2 - \a\, 
elsewise. 



(A7) 



The marginal distributions of a, G, and — a/G now 
follow by integration over P(a, G), 

P(a) = 1 - H, |a| < 1, (A8) 

P G (G) = 1 - \G- 1|, < G < 2, (A9) 

P(/3) = (l + |/3|)- 2 , \(3\<1, (A10) 



in accord with Eq. (18 1. We have plotted these distribu- 
tions in Fig. [3] 



0.5 



C 



t 



-1 1 

contact parameter 

FIG. 3: Probability distributions of the parameters at and 
/3i — cii/Gi that characterize a single-mode contact in the 
CRE, given by Eqs. |A8} and <A10| |. The distribution ( |A9 1 
of Gi — 1 is the same as that of a,:, but these two quantities 



(Al) are not independent because of the inequality &A6\ 



APPENDIX B: PROOF THAT THE TUNNEL 
RESISTANCE DROPS OUT OF THE NONLOCAL 
RESISTANCE 



According to Eq. ( 17 1, the nonlocal conductance G12 is 



determined by the product of the injection efficiency ol\ of 
contact N\ and the detection efficiency ct2 of contact N2. 
A tunnel barrier between the metal electrode and the su- 
perconductor suppresses the injection/detection efficien- 
cies and thereby suppresses the nonlocal conductance. 

The nonlocal resistance R12 is determined by the ra- 
tio oii/Gi of the injection/detection efficiency and the 
contact conductance Gi. Since both and Gi are sup- 
pressed by a tunnel barrier, one might hope that R12 
would remain of order e 2 /h. In this Appendix we in- 
vestigate the effect of a tunnel barrier on the nonlocal 
resistance, and demonstrate that it drops out identically 
for a single-mode contact between the normal metal and 
the superconductor. 

The key identity that we will use to prove this cancel- 
lation, is the product rule (fl6| and two corollaries: 



|Tr (Q Mi)a v (Q a y = \{{\Tv M % a v M?Oy) , 

i i i 

(Bla) 

\Tx{Ma y M T (jy)- 1 = [\TxMa y M T a y \~ 1 , (Bib) 

valid for arbitrary 2x2 matrices Mi. 

Considering any one of the two contacts, we assume 
that its scattering matrix So is modified by a tunnel bar- 
rier with scattering matrix 5S. Transmission and reflec- 
tion submatrices are defined as in Eq. (fsl) , 



S 



ro t 

¥ r' 
l '0 



6S 



Sr St 
5? Sr' 



(B2) 



For a single-mode contact, each submatrix has dimension 
2x2. Both So and SS are real orthogonal matrices at 
zero energy (in the basis of Majorana fermions). The 
tunnel barrier does not couple electrons and holes, which 
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means that the submatrices of 8S must commute with 

CTy, 



[cr y ,8r] = [(T y ,Sr'] = [a y ,8t] = [cr y ,8t'] = 



(B3) 



The submatrices of So are not so constrained. 

The total scattering matrix S of the contact is con- 
structed from 5*0 and 8S, according to the composition 
rule for scattering matrices. The transmission and reflec- 
tion submatrices of S take the form 

t = t (l -Srr'^St, (B4a) 

t' = St'(l - r' 5r)-h' , (B4b) 

r' = 8r' + 5t'r' (l - Srr'^St, (B4c) 

r = r Q + t 5r(l - r^r) -1 ^. (B4d) 

The injection efficiency a and detection efficiency a' 
are defined by 



^Trta y t T (jy, a' = |Tr t'a y t' T a y . 



(B5) 



Using the identities (Blal and (Bib I we can factor these 
quantities, 



a = ota5a/X, a = a' 6a' /X, 



(B6) 



into the product of the injection/detection efficiencies 
ao,a' without the tunnel barrier and terms containing 
the effect of the tunnel barrier: 

a = ^Tr t a y tQ(T y , a' — ^Tr t' a y t'Qcr y , (B7a) 
5a = sTr 8ta y 8t T cr y , 8a' = |Tr St' a y St' T a yi (B7b) 
X = |Tr (1 - 8rr' )a y (l - 5rr' ) T a y . (B7c) 
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simplify to 



Since 8t and 8t' commute with a y , the terms 8a, 8a' 



8a = 8a = ^Tr 8t8t T , 



(B8) 



where we have used the orthogonality condition, 
SS T SS = 8SSS T — 1, to equate the traces of 8tSt T and 
St' St' 1 . The term X can similarly be reduced to 

X = 1 + (1 G ) -Tr Srr' , (B9) 



where Go is the contact conductance (in units of e 2 /h) 
in the absence of the tunnel barrier: 



G = |Tr(l -r' a y r' a y ). 



(BIO) 



We now turn to the contact conductances G 2 ; , in order 
to show that the effect of the tunnel barrier is contained 
in the same factor Sa/X (which will then cancel out of 
the ratio f3i — ai/Gi). Considering again a single con- 
tact, and dropping the index i for ease of notation, we 
start from the definition of the contact conductance (in 
units of e 2 /h): 



G =JrTr (1 - r'a y r' T a y ). (Bll) 
We substitute Eq. ( B4c ) , and try to factor out the terms 



containing the transmission and reflection matrices of the 
tunnel barrier. 



It is helpful to first combine the two terms in Eq. ( B4c 
into a single term, using the orthogonality of SS: 



-{5t lT )- l 8r T 8t + St'r' (l - Srr'^-Ht 
(8t' T )-\r' - 8r T )(l - Srr'^St. 



(B12) 



We now substitute Eq. (B12) into Eq. (13 1 and use the 



identities ( B 1 ) to factor the trace 



G = 1 - X- 1 |Tr (r' - 8r T )<r y {r' T - 5r)a y 

= 1- X- 1 {2-5a-G -Tv5rr' Q ), (B13) 



where we also used the commutation relations (B3). The 



remaining trace of Srr' can be eliminated with the help 
of Eq. (B9l, and so we finally arrive at the desired result: 



G = G Q 8a/X. 



(B14) 



